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ABSTRACT 

The metallicity of a star strongly effects both its evolution and the properties of the stellar remnant that 
results from its demise. It is generally accepted that stars with initial masses below ~ 8 M© leave behind 
white dwarfs and that some sub-population of these lead to Type la supernovae. However, it is often tacitly 
assumed that metallicity has no effect on the rate of SNe la. We propose that a natural consequence of the 
effects of metallicity is to significantly increase the SN la rate in lower-metallicity galaxies. This is because 
lower-metallicity stars leave behind higher-mass white dwarfs, which should generally be easier to bring to an 
explosion. Using a simple model to relate the SN rate to galaxy age and metallicity, we find that the elevation 
in the rate of SNe la in lower-mass galaxies measured by LOSS is readily explained. We also find that models 
using the same parameters agree well with cosmic SN la rates up to z ps 2. We discuss additional implications 
of metallicity, including for inferences of the SN la delay time distribution and super-Chandrasekhar SNe. 
Subject headings: binaries: close — white dwarfs — supernovae: general 
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1. INTRODUCTION 

The end result of the evolution of stars that produce 
white dwarfs is often a Type la supernova explosion, 
through a single-degenerate channel (e.g., Whelan & Iben 
1973), double-degenerate channel (e.g., Iben & Tutukov 
1984; Webbink 1984), or both. Since stellar evolution is ob- 
viously affected by metallicity, there is no a priori reason why 
the rate of SNe la should not significantly depend on metal- 
licity. From a theoretical standpoint, a preference for high 
metallicity was proposed by Kobayashi et al. (1998), whose 
single-degenerate model required a minimum metallicity of 
<~ 0.1 Zq in order to produce Type la SNe. A similar pref- 
erence for higher metallicity was seen in the single degener- 
ate models of Langer et al. (2000). However, the strong pre- 
dictions offered by these models, such as no Type la SNe in 
dwarf galaxies and the outskirts of spirals, were not confirmed 
observationally (e.g., Prieto et al. 2008). The viability of the 
single-degenerate channel to produce the majority of Type la 
SNe has been debated from both observational (e.g., Leonard 
2007; Simon et al. 2009; Gilfanov & Bogdan 2010) and the- 
oretical viewpoints (Ruiter et al. 2009; Kasen et al. 2009; 
Hachisu et al. 2010), while simulations of double-degenerate 
mergers still do not yield normal SNe la (Pakmor et al. 2010; 
Fryer et al. 2010). 

From the viewpoint of stellar evolution, we expect an 
opposite sign for the dependence of the rate of Type la 
SNe on metallicity. Stars of lower metallicity at a given 
mass generally produce more massive white dwarfs accord- 
ing to stellar evolution calculations (e.g., Umeda et al. 1999; 
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Marigo & Girardi 2007; Meng et al. 2008), which should 
make it easier for them to reach the Chandrasekhar mass and 
explode. For example, according to Marigo & Girardi (2007), 
a single star with an initial mass of 3 Mq will leave behind 
a ~ 0.7 Mq white dwarf when evolved at solar metallicity 
(Zq = 0.019), while a star with the same initial mass at much 
lower metallicity (Z = 0.001) will leave behind a > 0.8 M Q 
white dwarf. Due to the steepness of the stellar initial mass 
function (IMF), this leads to a larger number of possible SN la 
progenitors. Obviously, producing a Type la SN explosion is 
a more complicated process than just evolving single stars (for 
example, see the discussion of common envelope phase treat- 
ment in Ruiter et al. 2009). However, the observed rate of 
Type la SNe implies that a large fraction (~ 2 — 40%) of all 
3 < M < 8 Mq stars will explode as one (e.g., Maoz 2008), 
which suggests that the binary evolution leading to SN la pro- 
duction cannot be "fragile". 

In this paper, we propose that the Type la supernova rate has 
a strong dependence on stellar metallicity and examine the re- 
sulting observational implications. This is motivated by mea- 
surements of the nearby Type la supernova rates reported by 
the Lick Observatory SN Search (LOSS: Leaman et al. 201 1; 
Li et al. 2011a,b). These measurements and their interpre- 
tation are discussed in Section 2. In Section 3, we present 
a simple model that examines the SN la rate as a function 
of galaxy mass/age, and incorporate the effects of metallic- 
ity in Section 4. In Section 5, we extend this into a treat- 
ment of the cosmic SN la rate, including the relative effects 
on "prompt" and "delayed" SNe. We summarize our argu- 
ments and include further discussion on analogous observa- 
tions, SNe la in galactic halos, galactic chemical evolution, 
and super-Chandrasekhar SNe in Section 6. 

2. THE TYPE IA SUPERNOVA RATE IN GALAXIES 

At present, the most complete and systematic search for 
nearby supernovae was conducted over the past decade by the 
Lick Observatory SN Search, with results recently detailed in 
Leaman et al. (201 1) and Li et al. (201 la,b). Here, we briefly 
discuss the implications of the LOSS findings for our present 
study. Of particular interest are the results pertaining to Type 
la supernovae. 

In Fig. 1, we display the specific SN la rate (rate per unit 
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FIG. 1 . — The specific rate of Type la supernovae versus host galaxy mass. 
Shown are data from LOSS for galaxies grouped by Hubble type (Li et al. 
2011b). Our models are also displayed, which assume either a At -1 de- 
lay time distribution alone (solid line) or an additional dependence on stellar 
metallicity (dashed, dotted lines; see text). 

mass) versus galaxy mass as measured by LOSS (Li et al. 
2011b). One is first struck by the steep dependence of this 
specific rate on galaxy mass. This variation of over an or- 
der of magnitude demands a physical explanation. The cause 
should be distinct from the origin of a similar pattern seen 
in the specific core-collapse supernova rate by LOSS, which 
likely arises mainly from the dependence of the specific star 
formation rate on galaxy mass (Li et al. 201 lb). 

Importantly, we also see that at a fixed mass the SN la rate 
does not vary greatly between galaxies of different Hubble 
type. This suggests that by examining a large set of galaxies 
one can arrive at the global behavior of SNe la. We proceed by 
translating the LOSS measured specific SN la rates in galax- 
ies of various Hubble types from a function of galactic mass 
into one of galactic metallicity. To do this, we convert be- 
tween galactic mass and median metallicity using the relation 
derived from SDSS data in Gallazzi et al. (2005), as shown by 
the upper axis of Fig. 1. This technique effectively averages 
over a large representative galaxy population similar to that 
sampled by LOSS. 

Since there is a delay from stellar birth to SN la explosion, 
a galaxy's SN la rate depends upon the age of its white dwarf 
population. This is typically quantified by an empirical or 
theoretical delay-time distribution (DTD), which results in a 
SN la rate that can be simply written as 



iV Ia (t) 



dt' 4>{t - t') p*0') 



(1) 



In 



where to is the age of the universe when SN la progenitor 
stars first formed and tp(t — t') is the DTD, which maps be- 
tween the rate of star formation at time t 1 , p*(t'), and the 
SN la rate at a later time t = t' + At. Eq. (1) can be used 
to calculate the expected SN la rate of an individual galaxy or 
the universe as a whole, given a properly normalized p* (t) 
(for the cosmic SN la rate, we will use the star formation 
rate density). Recent studies have suggested that <f> roughly 
takes a At -1 = (t - t'y 1 form (e.g., Totani et al. 2008; 



Maoz et al. 2010, 201 1). The physics behind this relation re- 
mains unclear, although such a distribution may naturally re- 
sult from binary mergers (see, e.g., Ruiteretal. 2009) or a 
single-degenerate scenario (Hachisu et al. 2008). 

Gallazzi et al. (2005) also derive r-band light-weighted 
galaxy ages, which vary from ~ 10 9 yr at 10 9 M Q to ~ 
10 10 yr at 10 12 M Q (see their Fig. 8), using galactic models 
with an exponentially declining star formation history from a 
time if orm with subsequent random bursts. Ideally, one would 
have at hand the detailed history of star formation in every 
galaxy. This is understandably difficult to achieve with any 
certainty. Attempts have been made in this direction (e.g., 
Brandt et al. 2010; Maoz et al. 2011); however, using what 
amounts to an average over the galaxy population should be 
suitable for comparison with global rates. We will use both of 
the Gallazzi et al. (2005) relations in estimating the expected 
SN la rate in galaxies of varying mass. 

3. A SIMPLE GALACTIC RATE MODEL 

We first attempt to explain the rate variations in galax- 
ies of different mass as due to an age effect alone. If the 
Gallazzi et al. (2005) ages corresponded to a single-age stel- 
lar population at a given galactic mass, deriving the expected 
SN la rate for a given DTD would be rather straightforward. 
For example, using a DTD for each galaxy of the form 



<t>{At) = 4>, 



(2) 



with tcyr = t/{\ Gyr), and assuming that the entire galactic 
stellar mass, M g , arose at a single time, t g , in Eq. (1) would 
lead to a galactic specific SN la rate at time t of 



N g (t) 
M„ 



<^Gyr - tg, Gyr) 



(t Gyr - t ffi Qyr) 7 



(3) 



This description is incomplete, though. First, the SN la rate 
at present reflects the galactic mass at the time of forma- 
tion, as opposed to that measured today after stellar mass 
loss has occurred. We correct for this using the results of 
Bruzual & Chariot (2003) for a Chabrier IMF (as used in 
the SDSS galaxy sample) by including a term of the form 
M(t g )/M(t). 

Additionally, the ages are more accurately galactic aver- 
ages, so that an assumption of instantaneous formation at t g 
will not properly reflect the effect of a DTD. To allow for a 
finite duration of star formation, we use a declining history of 
the form e _t / T , with r = 1 Gyr, occurring since the time t g 
for each galaxy. We further make use of the 16/84% ranges in 
log t g reported in Gallazzi et al. (2005) in order to weight the 
galaxy population with the DTD at fixed mass (rather than us- 
ing only the median value). These should alleviate the effect 
of average ages by giving greater weight to the low-age tail of 
their derived distribution and by allowing for a non-negligible 
rate of star formation today, particularly for lower-mass galax- 
ies, that is in rough agreement with the specific star formation 
rates measured by Schiminovich et al. (2007). 

This leaves the issue of the efficiency of converting a stel- 
lar population into SNe la (see, e.g., Maoz 2008). As we 
will discuss in Section 4, the DTD is also involved in shap- 
ing the cosmic SN la history, with again setting the overall 
normalization. Rather than attempting to incorporate theo- 
retical models of the DTD (see, e.g., Greggio 2010), we use 
a DTD with a pure power law of the form At -7 , with a 
lower cutoff t c to account for the minimum amount of time 
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needed to produce CO WDs. Evidence for delay times as 
short as < 100 Myr has been reported from, e.g., the study of 
SN remnants in the Magellanic Clouds (Badenes et al. 2009; 
Maoz & Badenes 2010), and we simply use t c = 50 Myr (see 
also the discussion in Section 5). 

The results of using this approach are shown as the solid 
line in Fig. 1, where we have used a At^^ DTD with 0* = 

1.4 x 10~ 3 (10 10 M ) _1 yr" 1 . A parametrization for this 
model is given in the Appendix. We see that the saturation 
in age at high masses results in a plateau, which should be a 
rather robust feature due to the relatively small scatter in esti- 
mated ages around < 10 Gyr in this range, while the decrease 
in age at lower mass results in a rise in the SN la rate. Overall, 
this simplified model agrees rather well with the LOSS data. 

4. INCORPORATING METALLICITY DEPENDENCE 

Historically, studies have focussed upon deriving the DTD 
without taking into account the possible effects of stellar 
metallicity on the SN la rate in a galaxy, as we have done 
above. As we next show, this assumption may not be suffi- 
cient to determine the actual DTD and to extract its astrophys- 
ical origins. We distinguish here between metallicity effects 
as primary (those involved in the rate of explosions) and sec- 
ondary (those affecting the detailed properties of individual 
explosions; e.g., Timmes et al. 2003), with our interest being 
in the former. 

We now examine a plausible scenario for including a SN la 
rate that varies with stellar metallicity. We propose that this 
arises from the effect of metallicity on the white dwarf pro- 
duced. In general, it is expected that, for the same initial 
stellar mass, the white dwarf from a star of lower metallic- 
ity should be more massive. This may be due to decreased 
mass loss and/or opacity resulting in hotter burning over the 
lifetime of the star (e.g., Umedaetal. 1999; Willson 2000; 
Marigo & Girardi 2007; Meng et al. 2008). The simplest in- 
terpretation of this is that it should be easier to reach the req- 
uisite Chandrasekhar mass for explosion through the addition 
of mass via binary evolution or a double white dwarf merger. 

To obtain a semi-quantitative estimate of the resulting 
change in the SN la rate with metallicity, we must consider 
the effect of a varying WD mass over the range of metallicities 
for the galaxies in the LOSS sample. One may hope for guid- 
ance from the initial-final WD mass relation determined from 
young star clusters. However, the clusters in which detailed 
studies are possible are nearby and formed recently, which 
necessarily limits them to single, ~ solar metallicity stars 
(e.g., Kalirai et al. 2008; Williams et al. 2009). We utilize the 
theoretical results of Umeda et al. (1999) in combination with 
the empirical metallicity estimates in Gallazzi et al. (2005). 
In decreasing the initial stellar metallicity from Z = 0.03 
to 0.004, roughly the range spanned in Fig. 1, Umedaetal. 
(1999) determined that an additional ~ 0.05 - 0.15 M© is 
added to the CO remnant at fixed initial mass (see their Fig. 6). 

Fig. 8 of Umedaetal. (1999) displays the relative num- 
ber of SN la progenitors obtained from their stellar evolu- 
tion model as a function of metallicity after integrating over 
a Salpeter IMF from a lower initial stellar mass (which varies 
with Z) corresponding to fixed final WD mass to an upper 
mass at which point ONeMg WDs were expected to be pro- 
duced (varying with Z from ~ 7—8.5 Mq). Using a threshold 
WD mass of 0.85 M Q yields a dependence on the rate with 
metallicity that can be approximately parametrized as 



over the range Z = 0.004 — 0.03. Using a lower threshold 
mass of 0.7 Mq yields a slightly weaker dependence, due to 
the larger mass range, of approximately 



iVi a (Z) oc (Z + 0.0015)~ 



(5) 



To derive galactic rates, we again use a SN la rate for each 
galaxy oc At" 7 and scale directly to the Umeda et al. (1999) 
results, normalizing these relations to unity at Z = 0.025, 
the metallicity of a characteristic ~ 10 11 Mq galaxy in the 
LOSS sample. Assuming Z g and t g to be separable, we use 
the 16/84% ranges in log Z g from Gallazzi et al. (2005) to 
again weight the galaxy distribution at fixed mass and intro- 
duce an overall term to account for the effect of metallicity in 
Eq. (1), either fom(Z) or f . 70 ( z )- 

The specific SN la rates resulting from using the two metal- 
licity scalings are shown in Fig. 1. For the Jo. 85 model, 0* = 



1.1 x 10- 3 (10 10 M Q )- ] 
while the /0.70 model has i 



yr" 1 with At G °; 8 (dashed line), 
>* = l.SxlO^ClO^A/o^yr- 1 
and At^y r 9 (dotted line). After accounting for the weaker ef- 
fects of the DTDs used, the rate does indeed rise more steeply 
at lower galactic masses than by taking into account age alone. 
This can be interpreted as a relative change in efficiency, an 
effect at the factor of ~ 2 level over the mass range of Fig. 1 
for the Jo. 85 model and slightly less for the /0.70 case. 

Because galaxy mass is strongly correlated with both age 
and metallicity, it is inevitable that the models are relatively 
degenerate and that inferences about the DTD from galaxy 
populations may err without accounting for metallicity. Our 
simplified treatment of galactic star formation histories may 
somewhat underpredict SN la rates at intermediate masses. 
This may be refined through more detailed modeling, al- 
though, given the uncertainties in our inputs, we will not at- 
tempt to do so here. It is encouraging that such broad agree- 
ment with data is already seen using quite general assump- 
tions. 

We note here that the normalization of these models can 
be scaled up or down, although this will directly affect the 
normalization of the expected cosmic SN la rates through the 
1. 
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FIG. 2. — The total masses of known WD-WD systems along with their 
calculated merger times due to gravitational wave losses. Shown are the three 
systems from the collection in Nelemans et al. (2005) in which both WDs 
have masses exceeding 0.5 Mq : two with firm masses (circles) and one with 
only a lower limit on the mass of the secondary (triangles). Assuming these 
to have all resulted from stars with Z ~ 0.02, we show the "expected" total 
masses for a range of metallicities (as labeled) using the final masses derived 
in Umedaetal. (1999). The resulting merger times assume initial orbital 
separations as presently inferred for each. 
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DTD, as we will discuss in the following Section. Note also 
that care should be taken in comparing these results, which 
examine the galaxy population as a whole, to those that dis- 
tinguish between "passive" and "star-forming" galaxies (e.g., 
Sullivan et al. 2006). 

The above is essentially based on an assumption of a single- 
degenerate scenario. While the full effects of metallicity on 
a double-degenerate scenario are likely more complicated, if 
we consider binaries with a uniform mass ratio distribution, 
the number that will exceed the Chandrasekhar limit depends 
on metallicity as 



N la (xZ 



O-l) b/a 



-0.4 



(6) 



where x — 2.35 is the slope of the IMF, a = 0.5 approxi- 
mates the slope of the white dwarf initial-final mass relation of 
Kalirai et al. (2008), and b ~ —0.08 is the dependence of the 
final mass on metallicity, estimated from Umeda et al. (1999). 
The magnitude of the effect is very similar to the case already 
considered, so we do not repeat the calculations. This model 
does not include any effect of the higher implied masses on 
the rate of binary evolution or possible effects in triple sys- 
tems based on the Kozai mechanism (Thompson 2010). 

To illustrate the above effect, we begin with the three known 
double WD binaries in which each component has a mass of 
at least 0.5 Mq from Nelemans et al. (2005). These have pri- 
mary/secondary masses of 0.71 + 0.55 Mq, 0.58 + 0.58 Mq, 
and 0.51 + >0.59 Mq, the last being a single-lined system 
with only a lower limit for the secondary. We assume that 
each WD arose from a solar-metallicity star and map from 
the WD masses to the initial stellar masses using the results 
of Umeda et al. (1999), mapping then to the WD masses cal- 
culated for these stellar masses at other metallicities. The 
0.5 Mq cut allows a straightforward translation without re- 
gards to systems with low-mass He WDs, etc. Keeping the 
initial binary separations fixed, we calculate the merger time 
due to gravitational wave losses for each system (see, e.g., 
Thompson 2010). Fig. 2 displays the effect on these systems 
using this prescription, where it is seen that the total masses of 
all three systems would have been pushed beyond the Chan- 
drasekhar mass limit and the merger time would have been 
significantly reduced at lower metallicities. 

There is hope for new tests to reduce the uncertainty in the 
overall effect of metallicity. For example, in observations of 
SN la host galaxies we would expect the hosts of SN la to be 
slightly less metal-rich than the galaxy population as a whole 
for fixed galaxy mass. This effect would not be as marked as 
in the case of gamma-ray bursts (see Stanek et al. 2006), since 
no hard upper metallicity threshold prohibiting the production 
of a SN la progenitor system is known to exist. An exception 
to this may be found at very high metallicity, as evidenced 
by the abundance of He rather than CO WDs in the metal- 
rich cluster NGC 6791 (see Kilic et al. 2007) - very massive, 
metal-rich galaxies may show an additional deficit of SNe la 
beyond that of our simple model. This may even be evident in 
the data at the high-mass end of Fig. 1, although it is difficult 
to draw a strong conclusion at present. 

We also note that the model that we have used only re- 
sults in a rather modest rate change with metallicity. It 
does not attempt to account for changes in the remnant mass 
that occur during the AGB phase of an isolated star (see, 
e.g., Vassiliadis & Wood 1993; Bird & Pinsonneault 2011; 
Renedo et al. 2010), which could result in a larger metallic- 
ity effect. Since both SN la scenarios require binary evolution 
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FIG. 3. — The cosmic rate of Type la supernovae. Shown are re- 
cent measurements from LOSS (Li et al. 2011b), SDSS (Dilday et al. 
2010), SCP (Kuznetsova et al. 2008), HST (Dahlen et al. 2008), SNLS 
(Gonzalez-Gaitan et al. 2011) and Subaru Deep Field (Graur et al. 2011). 
A model assuming only a fixed At^ 1 delay time distribution and the cos- 
mic SFR from Kistler et al. (2009) (thick solid line) can be compared to our 
models incorporating metallicity dependence (see text), which use either a 
At -0 ' 9 DTD (thick dotted line) or At' 8 DTD (thick dashed line). The 
components of these models with delays from stellar birth to explosion of 
less than 1 Gyr ("prompt") and greater than 1 Gyr ("delayed") are also shown 
(thin lines; as labelled). 



at some step in their evolution, this model should be adequate 
in this respect. We also have not attempted to vary the bi- 
nary fraction with stellar mass or metallicity. This is not yet 
well understood either theoretically or empirically (see, e.g., 
Mazeh et al. 2006), particularly in the mass range of the pro- 
genitors of SNe la, suggesting that substantial room for im- 
provement exists on both fronts in this regime. 

5. THE COSMIC TYPE IA SUPERNOVA RATE 

We next examine the expectations for the cosmic rate of 
SNe la by first considering a case without explicit metallicity 
dependence. We proceed by returning to Eq. (1) with the co- 
moving star formation rate density (z) inferred up to z ~ 8, 
using the Yuksel et al. (2008) parametrization of the star for- 
mation history (SFH), 



p*(z)=p 



(1 + z) Q " + 



1 + z 
B 



brj 



i 



c 



CI] 



1/'/ 

(7) 



where a = 3.4, b = —0.3, and c = —2, with breaks at Z\ = 1 
and Z2 = 4 corresponding to B = (1 + zi) 1-0 / 6 ~ 5100 
and C = (1 + z x )^- a ^ c {l + z 2 ) 1 " fc/c ~ 14, which re- 
flect the updated high-z data from Kistler et al. (2009), and 
we use r\ ~ —10 to smooth the transitions. The normaliza- 
tion is po = 0.014 Mq yr~ x Mpc~ 3 , which we have scaled 
down by a factor of 0.7 from the Salpeter IMF normaliza- 
tion of Hopkins & Beacom (2006) to be in better agreement 
with the galactic mass estimates in the previous Sections. In 
Eq. (1), we set to ~ 0.4 Gyr corresponding to z w 10. 9 

9 We use f2 m = 0.3, f!A = 0.7, and Ho = 70 km/s/Mpc where needed, 
e.g., in converting z o t and in rescaling the data in Fig. 3 using a common 
value of Hq. 
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Since we must consider the shortest possible delay times in 
constructing the cosmic rate from the SFH, we again assume 
a power law DTD with a cutoff at t c = 50 Myr. We note 
that the results of Umeda et al. (1999) (see also Siess 2007; 
Meng et al. 2008) suggest a maximum CO WD mass of ~ 
1.1 Mq that is nearly independent of metallicity. Since the 
effect of decreasing the metallicity is similar to increasing the 
stellar mass, we take this cutoff to be independent of Z (and 
thus z) since the lifetimes of the stars giving the most massive 
CO WDs should be similar. 

We use the same 0* and At^. DTD as in Section 3, with 
the resulting SN la rate history shown in Fig. 3 (thick solid 
line). To compare with prior results (e.g., Mannucci et al. 
2006; Scannapieco & Bildsten 2006; Sullivan et al. 2006), 
this history is broken down into the components with delay 
less than 1 Gyr (thin solid lines labelled as "prompt") and 
greater than 1 Gyr (labelled as "delayed"). We see that the 
"prompt" component is subdominant at z = 0, in agreement 
with the rates in Fig. 1. Altering either the form of the DTD or 
t c can make the "prompt" component relatively more or less 
important (see, e.g., Horiuchi & Beacom 2010 for related dis- 
cussion); however, this would in turn affect the specific SN la 
rate models in Sections 3 and 4. 

Since the universe as whole had a lower metallicity in the 
past, a relative enhancement should also be effected in the 
cosmic SN la rate. As existing rate measurements average 
over the entirety of the galaxy population, this effect should 
not be dramatic at the present epoch, but, as for the spe- 
cific rate, can be important in deriving the DTD. At low 
z, the gas-phase metallicity is typically higher than that of 
the stellar population (Gallazzi et al. 2005). The relation be- 
tween galaxy mass and gas-phase metallicity is well deter- 
mined at low z (Tremonti et al. 2004) and has been measured 
to evolve at higher redshifts, so that the typical metallicity 
decreases by ~ 0.15 dex per z up to at least z « 2 (e.g., 
Kewley & Kobulnicky 2005; Savaglio et al. 2005; Erbetal. 
2006; Maiolino et al. 2008). We use 

Z(z) = 0.03 x lCT al5z (8) 

to account for stars forming from gas that is increasingly 
metal poor at higher z, with a resulting change in rate aris- 
ing through either the relation approximated by Eq. (4) for 
the no. 85 model or Eq. (5) for the ?io.70 model, again normal- 
izing each to unity at Z = 0.025 to be consistent with our 
specific rate models. We also use the same values of (f>* and 
DTD slopes as in the corresponding specific rate models. 

Fig. 3 shows the resulting cosmic rates for the tt-o.85 (with 
AiQ° r 8 ; thick dashed line) and n .7o (Aic° r 9 ; thick dot- 
ted line) models. Both models yield similar histories as the 
metallicity-independent case, with parametrizations for all 
three included in the Appendix. This is due to the relative 
increase of the rates with z as compared to models with the 
same DTD without a metallicity enhancement. This is similar 
in spirit, but less dramatic, than the relative evolution likely 
due to stellar metallicity seen in the cosmic GRB rate (e.g., 
Kistler et al. 2008). Both models are also broken down by de- 
lay time in Fig. 3 (thin dotted, dashed lines), which illustrates 
the underlying effect of altering the DTD. 

As discussed for the specific SN la rate, there is again a de- 
generacy between altering the DTD and including the effect 
of metallicity, although not quite as strong. That the metallic- 
ity effect works in the same direction as decreasing the index 
in the DTD in both cases, as seen in Fig. 1 and Fig. 3, is some- 



thing of a coincidence, owing to the fact that galaxy ages and 
metallicities both decrease with decreasing mass and the cos- 
mic SFR rises with increasing z. This didn't have to be the 
case, though. We thus reiterate that an estimate of one com- 
ponent must account for the other until this degeneracy is bro- 
ken. It is possible to perform a more elaborate study by vary- 
ing all the parameters involved (see, e.g., Horiuchi & Beacom 
2010 and Grauretal. 2011 for the metallicity-independent 
DTD); however, the qualitative effects of the models that we 
have considered are already sufficiently evident. 

6. DISCUSSION AND CONCLUSIONS 

The rate of Type la supernovae should be affected at some 
level by the effects of metallicity on stellar evolution. There 
may be various complications involved, such as the largely 
unresolved effects of binary evolution, but our simple model 
for the effects of metallicity should be broadly relevant. There 
has been significant effort devoted to investigating Type la 
SN properties as a function of metallicity (e.g., Hamuy et al. 
2000; Gallagher et al. 2008; Howell et al. 2009; Neill et al. 
2009; Sullivan et al. 2010; Konishi et al. 2011). Since the 
properties of SN la have been seen to depend on metallicity, 
why not the rate? 

The simple models that we have considered explain fairly 
well both the specific SN la rates measured in nearby galax- 
ies by LOSS and the observed normalization and evolution of 
the cosmic SN la rate. Further progress can certainly be made 
with data that can break the degeneracy between decreasing 
metallicity and decreasing age. Indeed, evidence in this di- 
rection has been found in a comparison of SN la host galaxies 
in SDSS by Cooper et al. (2009), who found that SNe la in 
blue, star-forming hosts have a preference for lower-density 
environments, which they interpreted as being the effect of 
lower gas-phase metallicities. 

As previously mentioned, the results of Umeda et al. (1999) 
indicate that the maximum CO WD mass remains close to 
~ 1.1 Mq over a wide range of metallicities. If this is true, 
and binary evolution effects are neglected, then we would ex- 
pect the relative rates of super-Chandrasekhar SNe la arising 
from mergers to increase with lower metallicity in proportion 
to the normal SNe la due to the power law form of the IMF 
(note that instabilities prohibit the necessary growth of even 
rapidly-rotating single WDs; Piro 2008). However, recent 
observations of host galaxies may indicate an even stronger 
preference for low-metallicity hosts for super-Chandrasekhar 
SNe la (e.g., Taubenberger et al. 2011; Childress et al. 2011). 
Moreover, the maximum mass resulting from a merger under 
these assumptions is ~ 2.2 Mq, well below that inferred from 
SN 2007if (Scalzo et al. 2010). 

An explanation for both effects may arise from binary evo- 
lution. To achieve a higher total merger mass, without re- 
sorting to an ONeMg WD, at least one WD should gain 
mass while maintaining a CO composition. If the effect 
of inhibiting single degenerate SN la production at Z < 
0.1 Z Q (Kobayashi et al. 1998; Kobayashi & Nomoto 2009) 
does hold, then the primary WD may instead be pushed close 
to, but not above, the threshold for explosion, so that the rate 
of massive mergers is further enhanced. The end state of the 
secondary resulting in a massive CO WD could then lead to a 
merger with total mass upwards of ~ 2.5 Mq. 

If, contrary to Umeda et al. (1999), the CO WD mass limit 
actually increases modestly for lower metallicities, then the 
rate of such extreme super-Chandrasekhar mergers rises dra- 
matically at lower metallicities without a need to turn to bi- 
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nary evolution for a solution. This is due to the presence of the 
threshold in reaching the requisite total merger mass, which 
would lead to a large relative difference between low/high- 
Z galaxies. The stars giving rise to these massive white 
dwarfs would also evolve more rapidly and could thus lead 
to "prompter" explosions in low-Z environments. Whether 
this scenario occurs is a question left for stellar evolutionary 
modeling and observations of host galaxies. 

Additionally, observations of Milky Way dwarf spheroidal 
galaxies have revealed decreasing values of [a/Fe] with in- 
creasing [Fe/H], indicating the influence of Type la super- 
novae down to metallicities of [Fe/H] « —2.5 (Kirby et al. 
2011). Since model fits to these measurements are naturally 
sensitive to the SN la rate over a range of metallicities, we 
urge exploration of the implications of an increased rate at 
low Z, including super-Chandrasekhar mergers, on galactic 
chemical evolution. 

In the category of interesting, but more anecdotal, evidence 
that low metallicity might be of significance for Type la su- 
pernovae, Tovmassian et al. (2010) recently presented strong 
evidence that SBS1 150+599 A, a close binary star inside a 
metal-poor, Galactic halo planetary nebula PN G135.9+55.9 
consists of two white dwarfs that will merge within a Hubble 
time. The estimated total mass of the binary is very close to 
the Chandrasekhar limit, making it a likely SN la progenitor. 

It is also interesting to note that the normalized rate of 
planetary nebulae in elliptical galaxies (Buzzoni et al. 2006) 



shows a very similar trend with metallicity to that discussed 
in our Fig. 1. Indeed, their Figs. 11 and 12 show about 10 
times fewer planetary nebulae per unit luminosity in metal- 
rich, massive ellipticals compared to metal-poor, low-mass 
ellipticals. The mapping between PN production and SN la 
explosion is of course uncertain; however, both involve the 
production of a white dwarf, and Buzzoni et al. (2006) at- 
tribute finding fewer PNe in more metal rich ellipticals to a 
dependence of the initial-to-final mass relation on metallicity. 

Substantial observational progress has been made in the 
study of SNe la in the last decade and new data can be ex- 
pected to better determine the extent to which metallicity af- 
fects the SN la rate. Possible directions include detailed mea- 
surements of rates as a function of galaxy mass. Related to 
this is the expectation of a relatively higher SN la rate in 
galactic halos, for a given stellar population age and total 
mass, due to the lower average metallicity. Such a bias may 
already be seen in the number of cases in which a SN la oc- 
curred in the outer halo of a star-forming galaxy (Prieto et al. 
2008; Khan et al. 201 1). Elevated rates of other transients po- 
tentially involving a WD (e.g., Thompson et al. 2009) should 
also be expected. Improved measurements of the cosmic 
SN la rate, in combination with independent determinations 
of the DTD at constant metallicity, can determine whether the 
rate is larger than otherwise expected. Along with this, if the 
intrinsic properties of SNe la vary with metallicity, we would 
expect an evolution in the Type la luminosity function as the 
lower- Z component becomes prominent at higher redshifts. 



APPENDIX 



The models that we have discussed are the result of combining several unrelated functions and thus are not necessarily 
amenable to convenient parametrization. Nonetheless, we find that a sigmoid function provides an adequate fit to our metallicity - 
independent model of the specific SN la rate, (j a , with 



Cla(Af) 



(1O 1O M )-V- 



1 + cxp 



log(M/M Q ) - Ah 



(1) 



where q = 5x 10~ 3 , (3 = 4.2 x 10 -4 , M* = 10, and u> = 0.33 agrees with the model to within < 10% over the mass range 
displayed in Fig. 1. The metallicity-dependent models can be fit with similar parameters. 

Using the smoothly-broken piecewise form of Eq. (7), with po replaced by ho, our cosmic rate models can be fit to within a 
few percent over the range z = — 4. The parameters used for the metallicity-independent model and the metallicity-dependent 
no. 85 and no. 70 models are given below in Table 1. All three use r\ ~ —10 to smooth the transitions. 



TABLE 1 

Parameters used in the fits of our three cosmic SN Ia rate scenarios. 



Model 


no [yr 1 Mpc 3 ] 


a 


b 


c 


Zl 


Z2 


B 


C 


Z-free 


2.5 X 1CT 5 


1.8 


-0.8 


-2.3 


0.9 


2.9 


8.1 


5.0 


"0.85 


2.9 X 10" 5 


1.4 


-0.5 


-2.0 


0.9 


2.9 


11.5 


5.1 


"0.70 


2.8 x IO- 5 


1.6 


-0.7 


-2.0 


0.9 


3.0 


8.2 


5.2 
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